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the latest Kami AND results 
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Abstract. We present an updated analysis of all available solar and reactor neutrino data, emphasizing in 
particular the totality of the KamLAND (314d live time) results and including for the first time the solar 
SNO (391d live time, phase II NaCl-enhanced) spectrum data. As a novelty of the statistical analysis, we 
study the variability of the KamLand results with respect the use of diverse statistics. A new statistic, 
not used before is proposed. Moreover, in the analysis of the SNO spectrum a novel technique is used in 
order to include full correlated errors among bins. Combining all data, we obtain the following best-fit 
parameters: we determine individual neutrino mixing parameters and their errors Am^ — 8.2 ± 0.08 x 



10-^eV^ 



tan 
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elements of the neutrino mass matrix, 
matrix: 



The impact of these results is discussed. We also estimate the individual 
In the framework of three neutrino oscillations we obtain the mass 



M 



/ 1.0 + 4.0 ±3.2 10" 
eV 4.2 ±3.2 10"^ 
V 13.5 ± 14.0 10"^ 



4.2 ±3.2 10"" 
1.0 + 4.3 + 3.5 10"^ 
-13.5 + 14.5 10"^ 



-13.5 ± 14.0 lO"'' 
-13.5 ± 14.5 10"^ 
1.0 + 100.0 + 30.0 10 



PACS. 2 6.65.+t, 14.60.Pq 



1 Introduction 

Evidence of antineutrino disappearance in a beam of an- 
tineutrinos in the Kamland experiment has been presented 
[1]. The analysis of previous experimental results on reac- 
tor physics and solar neutrinos [2] in terms of neutrino 
oscillations has largely improved our knowledge of neu- 
trino mixing. Thus, the solar neutrino data evidence prior 
to autumn of 2003 converged to two distinct allowed re- 
gions in parameter space, often referred to as LMAI (cen- 
tered around the best-fit point of Aitiq = 7.1 x 10^''^eV^, 
tan^ 6*0 = 0.47) and LMAII (centered around Z\rn| = 
1.5 X 10-''eV^ tan^eQ = 0.48). The inclusion of SNO 
phase II data eliminated the LMAII region at about 4cr. 

The KamLAND measurementes presented recently [1] 
corresponds to the first about 317 days live time, they 
confirm the conclusions obtained from previous data and 
give a more stringent limit on the neutrino mass param- 
eters. As we will see in this work, the new KamLAND 
information further excludes the LMAII region now at sa 
5a. 

In addition to an increased statistics, a significant change 
in the KamLAND analysis technique is related to the 



fiducial volume definition. Whereas in the previous setup, 
events taking place at the outer edge of the nylon bal- 
loon were rejected [1,3], the recent analysis adopts a more 
sophisticated coincidence-measurement technique to ex- 
clude unwanted backgrounds. Additionally, a better un- 
derstanding of the fuel cycle on the reactors has lead to 
the collaboration to estimate the incoming neutrino flux 
with better accuracy: the estimated error on this initial 
flux (po (Ve) is now of the order of 2%. 

The aim of this work is to present a comprehensive up- 
dated analysis of all recent solar neutrino data including 
the KamLAND reactor-experiment results to determine 
the extent of the remaining viable region in the parameter 
space and to obtain favoured values for the neutrino phys- 
ical parameters in a two-neutrino framework, and, with 
the inclusion of results coming from atmospheric oscilla- 
tion evidence give an estimation of the elements of the 
three neutrino mass matrix. Some key analysis novelties 
are presented along this work, for example regarding the 
treatment of sistematic correlations in the analysis of the 
SNO spectrum (day+night) and a improved analysis of 
the KL data with an appropriate consideration of its low 
statistics data bins. 
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The structure of this paper is the following: in section 2 
discuss our approach to the latest KamLAND results and 
older solar evidence. All solar neutrino experiments are 
discussed in section 2.1. We specially discuss the impor- 
tance of the SNO data and the spectrum results. We then 
proceed, in section 3, to explain the procedure adopted 
in our analysis with some emphasys on the special treat- 
ment needed by the KamLand and SNO spectra and in 
section 4 we present our results. Finally we summarize 
and conclude in section 5. 



2 The KamLAND and solar evidence 

Reactor anti- neutrinos with energies above 1.8 Mc^V pro- 
duced in some 53 commercial reactors arc detected in 
the KamLAND detector via the inverse /3-decay reaction 
Vp +p^ n + e+ . The mean reactor-detector distance and 
energy window of these Vg makes KamLAND an ideal test- 
ing ground for the LMA region of the Vq parameter space. 
The first results published by the KamLAND collabora- 
tion eliminated all possible solutions to the solar neutrino 
problem (SNP) except the LMA region of the parameter 
space [3]. The sensitivity of this experiment to the Aw? 
parameter divided the previously whole LMA region into 
two distinct regions, the one relative to the smaller mass- 
squared difference being preferred by data [4]. 

In our analysis (see also Ref.[17] for further details) 
we model the reactor incoming flux through a constant, 
time-averaged fuel composition for all of the commercial 
reactors within detectable distance of the Kamioka site, 
namely "^^^U = 56.3%, 238[/ = 7.9%, asQp^ ^ 30.1%^ and 
^'^^Pu = 5.7%. We used the full cross-section including 
electron recoil corrections. We analyzed the data above 
threshold of 2.6MeV, as the low-energy end of the spec- 
trum had effectively no events. The information relative 
to the no-oscillation initial flux for the two periods can 
be extracted from fig. (l.a) of [1]. We neglect all back- 
grounds, including geological background above 2.6 MeV. 
We use the resolutions published by the collaboration for 
the two different data sets, namely a (E) = 6.2%/ y/E for 
the recent data (post upgrade) and cr [E] = 7.3%/a/B for 
earlier data (pre- upgrade) . The total systematic error is 
estimated at 6.5%. 

In order to use all the data available, we use a sim- 
ple MC simulation to estimate an equivalent efficiency for 
the two pre-upgrade and post-upgrade phases. Finally, no 
matter efii'ects are taken into consideration for the Kam- 
LAND data alone, as it was shown that, for this experi- 
ment, any asymmetry due to matter effects is negligible 
for Aw? of the order of 10 x 10^^ eV^. 



2.1 Solar data 

The most ponderous data present in our analysis come 
from the solar neutrino experiments. The experimental 
results are compared to an expected signal which we com- 
pute numerically by convoluting solar neutrino fluxes [5], 



sun and earth oscillation probabilities, neutrino cross sec- 
tions and detector energy response functions. We closely 
follow the same methods already well explained in pre- 
vious works [10, 12 14], we will mention here only a few 
aspects of this computation. We determine the neutrino 
oscillation probabilities using the standard methods found 
in literature [17], as explained in detail in [10] and in [12]. 
We use a thoroughly numerical method to calculate the 
neutrino evolution equations in the presence of matter 
for all the parameter space. For the solar neutrino case 
the calculation is split in three steps, corresponding to 
the neutrino propagation inside the Sun, in the vacuum 
(where the propagation is computed analytically) and in 
the Earth. We average over the neutrino production point 
inside the Sun and we take the electron number density 
rie in the Sun by the BPB2001 model [5]. The averaging 
over the annual variation of the orbit is also exactly per- 
formed. To take the Earth matter effects into account, we 
adopt a spherical model of the Earth density and chemical 
composition [17]. The joining of the neutrino propagation 
in the three different regions is performed exactly using 
an evolution operator formalism [17]. The final survival 
probabilities are obtained from the corresponding (non- 
pure) density matrices built from the evolution operators 
in each of these three regions. 

In summary, double-binned day-night and zenith an- 
gle bins are computed in order to analyze the full Su- 
perKamiokande data [7], whereas single-binned data is 
used for the SNO detector [8,16,18]. The global signals 
only are used for the radiochemical experiments Homes- 
take [19], SAGE [20,21], GallEx [22] and GNO [23]. The 
next paragraphs are dedicated to a description of SNO 
component of the solar evidence. 

The Sudbury Neutrino Observatory (SNO) collabora- 
tion has presented data relative to the NaCl phase of the 
experiment [8, 9). We remind that the addition of NaCl to 
a pure D2O detection medium has the effect of increas- 
ing the detector's sensitivity to the neutral-current (NC) 
reactions within its fiducial volume. The NC detection ef- 
ficiency has changed from a previous 'no-salt' phase of ap- 
proximately a factor three. This and other novelties have 
made it possible for the SNO collaboration to analyze their 
data without making use of the no-spectrum-distortion 
hypothesis. Furthermore, they have adopted a new 'event 
topology' criterion [9] to distinguish among the different 
channels within the detector. The SNO Collaboration has 
now devised a new data-analysis technique which relies 
on the topology of the three different events. The new pa- 
rameter {(3e) relative to which they marginalize is known 
as the 'isotropy' of the Cerenkov light distribution was 
used to separate the CC, ES and NC signals, something 
that was not possible in th previous two data sets. The 
measured fluxes are reported in [8]. 

The comparison of the new SNO results and the previ- 
ous phase-II data can easily be made because the SNO col- 
laboration has included in the recent paper results which 
were obtained following their previous method, along with 
the new unconstrained data. The new results are compat- 
ible with the previous ones. It seems that the overall ef- 
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feet of un-eonstraining the analysis is an increase in the 
measured fluxes, although the estimated total has de- 
creased relative to the previous best-flt value, leaving even 
less space for eventual sterile neutrinos. Wc make use of 
all day+night published data and refer the reader to sec- 
tion 3 for details. The procedure we used to introduce the 
SNO spectrum data is an extention of the one use for the 
SK spectrum analysis and is explained in the next section. 



3 Our analysis 

Wc use standard statistical techniques to test the non- 
oscillation hypothesis. Two different sets of analyses are 
possible with the present data on neutrino oscillations: 1) 
short-baseline reactor data, solar data including the SK 
spectrum and previous phase-I (CC only) SNO spectrum, 
phase-II SNO global result, combined with new the Kam- 
LAND spectrum and, 2) the previous set with the use 
of the phase-II SNO spectrum result and the new Kam- 
LAND data. In order to use all the SNO data, we consider 
the phase-I and phase-II results as two distinct but fully 
correlated experiments. 

For the purpose of this analysis, a function is de- 
fined which is the the sum of the distinct contributions. 
The contribution of all solar neutrino experiments is sum- 
marized in the term; 

Xlun = Xrad + XSK + XSNO- (1) 

Where, the function for the global rates of the radio- 
chemical experiments is as follows, 

(2) 

where R*"'exp g^j.g length-two vectors containing the the- 
oretical (or experimental) signal-to-no-oscillation expec- 
tation for the chlorine and gallium experiments. Corre- 
lated systematic and uncorrelated statistical errors are 
considered in agyst and astat respectively. Note that the 
parameter-dependent i?*'* is an averaged day-night quan- 
tity, as the radiochemical experiments are not sensitive to 
day-night variations. For the next component We consider 
the double-binned SK spectrum comprising of 8 energy 
bins for a total of 6 night bins and one day bin. The is 
given by 

xIk = ("R*" - R^^^)^ (^Lc + ^corV' ("R*'^ - R^^'')- 

(3) 

The covariance matrix tr is a 4-rank tensor containing 
information relative to the statistical errors and energy 
and zenith-angle bin-correlated and uncorrelated uncer- 
tainties. Since the publication of the first SNO NC results, 
we have adopted their estimate oi 4>b and incorporated 
the new parameter a in the representing the normal- 
ization with respect to this quantity. In determining our 
best-fit points, we minimize with respect to it. Note that 
the quantities R«^p and R*'' contain the number of events 
normalized to the no-oscillation scenario. 



We deal next with the SNO component. We present 
two different analysis of some of the phase II SNO data 
sets including total day/night quantities. 

The first consideres the global signal alone, the second 
incorporates the total signal spectrum. We consider here 
the two SNO results as if coming from two independent 
experiments, but fully correlated. We use the backgrounds 
as listed in tables X of Ref.[8] and 1 of Ref.[18] for the 
phase-II data. The detector resolution is obtained from 
Refs. [8,27]. 

In the second case we make use of the spectral data. 

The spectrum used for our analusis is presented in ta- 
ble [2]). The XsNO same formal expression as be- 
fore, where it is understood that R*'*'^^? are now length 32 
containing two 16-bin relative to the two SNO data sets. 
We consider the two fully correlated. 

The main difficulty in using the total spectrum data 
lies in correctly estimating the, highly correlated, system- 
atic error. By using the information contained in tables 
XIX and XX of Ref. [8] , we have computed the influence of 
all the different sources of error on our response function 
considering the correlation/anti-corrclation as presented 
in table 1 of [27]. The different backgrounds spectral cor- 
relations are included from table XXXIV of Ref. [8]. The 
procedure we used to introduce the SNO spectrum data 
is an extension of the one used for the SK spectrum anal- 
ysis. For each point in the parameter space Z\m^, 6 we 
start from a correlation matrix obtained by using the non- 
deformed spectrum assumption. We calculate, for each 
bin, the sum of the signals ES-|-NC-|-CC and we extract 
weights for each single contributions. After that, we com- 
pare our theoretical results with the ones given by the 
SNO collaboration and impose a 3-ct cut. By using these 
zero-order weights and the correlation errors obtained by 
the SNO table, we reconstruct a correlation matrix. The 
correlation matrix is introduced into the x^ analysis by 
adding a free parameter 5 cor which is determined in a min- 
imization process togethed with the weights of the single 
i = ES + NC + CC contributions to the signal: 

i 

The full process is designed to be interated a number 
of times, in practise we obtain that after two iterations 
the process is convergent and give us the desired results. 

3.1 The Kamland statistical analysis 

The total KamLAND contribution to the is deflned as: 

Xkl = Xkl ,glob + Xkl,x (4) 
where the global contribution is simply 

XKL,glob — ~2 • (^) 

^ stat+sys 

The statistical consideration of the KL spectrum sig- 
nal, xlfL A' is however worthy of special attention. Due to 
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the fact that at high energy KamLAND observes a small 
number of events alternatives should be used instead the 
Gaussian approximation. This means among other things 
that the correlated systematic deviations cannot be in- 
troduced in a straightforward way. Due to these reasons, 
we use an alternative technique for the KamLAND data 
bins. A detailed account of some statistical considerations 
is presented in the Appendix. 

It is possible to present an unified approach [29] to all 
the commonly used multinomial models (Pearson's, log- 
likelihood among them) by defining a family of statistics 
X^(A) for testing the fit of observed frequencies i?^^^ to the 
expected ones i?*'' [29] . All the statistics belonging to this 
family have similar well-behaved properties but however 
results as best fit parameters and exclusion regions may 
significantly depend on the use of one or another. Any 
decision as to which member of the family we should use 
to finally test the null hypothesis must depend on the type 
of the departure we wish to detect. The sensitivity of the 
statistic depends on how the defining function treats the 
large or small deviations. 

Based on a comparative study it is recomended [28, 29] 
to use x^(2/3) as a compromise candidate among the dif- 
ferent test statistics optimizing diverse criteria as rate of 
convergence, sensitivity to the sample size and sensitivity 
to large or small bin deviations. The statistic correspond- 
ing to this value, the Read statistic, is the one used in this 
work: 




In the evaluation of Xkl a '^^ vectors that comprise 
therefore of 13 spectral points of width 0.425 MeV. 



4 Results and Discussion. 

To test a particular oscillation hypothesis against the pa- 
rameters of the best fit (null hypothesis) and obtain al- 
lowed regions in parameter space we perform a minimisa- 
tion of the full function with respect the oscillation and 
the rest of ancillary parameters. A given point in the oscil- 
lation parameter space is allowed if the globally subtracted 
quantity fulfills the condition Ax^ = x^{Am?, ^) — Xmin < 
Xn{CL). Where X(90%,95%, ...) are the respective quan- 
tiles. In this way we obtain best fit mass differences and 
angles and joint exclusion regions. Additionally, we per- 
form a second kind of analysis in order to obtain concrete 
values for the individual oscillation parameters and esti- 
mates for their uncertainties. We study the marginalised 
parameter constraints where the x^ quantity is converted 
into likelihood using the expression £//!o = e~^^ -Xmi„)/2 
In table [4] we report the values of the mixing param- 
eters ArriQ, tan^ 9q, and the x^ obtained from minimiza- 
tion and from the peak of marginal likelihood distribution. 



The results are shown in Figs.l where we have gc!n- 
erated acceptance contours in the Avn?-t&v? 9 plane. In 
fig. [l-(left)] we show the exclusion plots for the solar, ra- 
diochemical + Cerenkov solar data and KamLAND with 
the global signal of the SNO phase-II data, whereas the 
right panel refers to the KamLAND spectrum, radiochem- 
ical -|- Cerenkov solar data and the SNO phase-II spectrum 
information. Contour lines correspond to the the allowed 
areas at 90, 95, 99 and 99.7% CL relative to the absolute 
minimum. 

Thes normalized marginal likelihood, obtained from 

the integration of C for each of the variables, is plotted in 
Figs. (2) for each of the oscillation parameters Am? and 
tan^ 9. Concrete values for the parameters are extracted 
by fitting one- or two-sided Gaussian distributions to any 
of the peaks (fits not showed in the plots). In both cases, 
for angle and the mass difference distributions the good- 
ness of fit of the Gaussian fit to each individual peak is 
excellent (g.o.f ~ 100%). The errors obtained from this 
method arc assigned to the x^ minimisation values. The 
central values are fully consistent and very similar to the 
values obtained from simple x^ minimisation. Systemat- 
ics variability of these results can come from the use of a 
different prior information or mixing parameterizations, 
however this variability or systematic error due to the 
procedure is small. We will again use the technique of 
marginal distributions in the next paragraphs to obtain 
an estimation of the individual elements of the neutrino 
mass matrix and their errors. 

The main difference with previous analysis is a better 
resolution in parameter space. The previously two well 
separated solutions LMAI,LMAII have now completely 
disappeared. In particular the secondary region at larger 
mass differences (LMAII) is now completely excluded. 

The introduction of the new KamLand data in general 
strongly diminishes the favored value for the mixing an- 
gle with respect to the KamLAND result alone [15]. The 
final value is more near to those values favored by the 
solar data alone than to the KamLAND ones. As an im- 
portant consequence, the combined analysis of solar and 
KamLAND data concludes that maximal mixing is not 
favored at ~ 4 — 5cr. This conclusion is not supported 
by the antineutrino, earth-controlled, conceptually sim- 
pler KamLAND results alone. As we already pointed out 
in Ref. [15], this effect could be simply due to the present 
low KamLAND statistics or, more worrying, to some sta- 
tistical artifact derived from the complexity of the analysis 
and of the heterogeneity of binned data involved. 

4.1 An estimation of the neutrino mass matrix 

We proceed now to an estimation of the neutrino mass 
matrix in different aproximations. Our main objective is 
however to estimate how well the individual errors of the 
mass matrix can be extracted already at present by the 
existing experimental evidence. For this purpose we have 
applied similar arguments as those used before to obtain 
marginal distributions and errors for individual parame- 
ters from them. 
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The square of the neutrino mass matrix can be written 
in the flavour basis as = UMj^W where Md is diag- 
onal and U is an unitary (purely active oscillations are 
assumed) mixing matrix. Subtracting one of the diagonal 
entries we have 

= m\l + Ml = mil + UM'^U\ 

where I is the identity matrix. In this way we distinguish 
in the mass matrix a part, A'/q, which affects and can be 
determined by oscillation experiments and another one, 
m\I, which does not. Evidently, the off-diagonal elements 
of the mass matrix are fully measurable by oscillation ex- 
periments. 

First, we restric ourselves for the sake of simplicity to 
two neutrino oscillations, we have in this case 



= mil + M,? = mil + Am^ 



sin^ 



sin 9 cos 6* \ 
,2 a (6) 



sin d cos d cos 



with Am? = m| — ml. The individual elements of the ma- 
trix Mo can simply be estimated from the oscillation pa- 
rameters obtained before. For example for tan^ 9 ~ 0.40, 
Am^ ~ 7 - 8 X 10-5 eV^ we would obtain (M^)22 ~ 
5 - 6 X 10-5 ey2. 

Using again as likelihood function the quantity 
£/£o(^rn2,tan2 6*) = e-^^"" -^^r^.J n ^e obtained the indi- 
vidual probability distributions for any of the elements of 
the matrix A^q. Average values and Icr errors arc obtained 
from two-sided Gaussian fits to these distributions. From 
this procedure we obtain: 



'2.05t°J5 

o 19-1-0.25 



o -,9+0.25 ^ 
"^•^^-0.26 

4.50t°iJ. 



(7) 



One can go further supposing a concrete value for m\ 
from elsewhere. If we take ml » Am^ then we can di- 
rectly write the mass matrix 



M = mil + 



1 



2m 1 



(8) 



Supposing for example mi 

2+o;2 10-5 



= 1 eV, 



M-.vf^-^ + 1-021^12 10-' i-56+°:ll 10-' \ 

^-^^y 1.56+° ll 10-' 1-0 + 2.25+°ll 10-' J ■ 

this is the final two neutrino mass matrix which can be 
obtained from present oscillation evidence coming from 
solar and reactor neutrinos. 

We obtain now an estimation of the three neutrino 
mass matrix. For this purpose we make the same reasoning 
as before and introduce the existing evidence of the indi- 
vidual values of the two additional angles and the square 
mass difference. Naturally knowledge of these parameters 
is still very poor and the elements of the final mass matrix 
will have much larger errors. First we substract a diagonal 
part and write the square mass matrix as: 



= mil ■ 



Aml^Ml^ 



Aml^Ml^ 



with Am^j = m^ — m'j. 



We write the mixing matrix as a product of three single 
rotations around each of the axis: 

U = Ul2(6'l2)W23(^'23)■Wl3(6'l3)• 
With this notation the matrices M12, M32 are written 

M12 = {U23U13YMQU23U13, 

M32 = iU23Ul3Y M3U23U13 

where M3 = Diag{0,0, 1). The matrix M32 does not de- 
pend on the angle ^12. The dependence on this angle is 
fully contained in Mq which is the 3x3 enlarged version 

of the Mq matrix appearing in Eqs.(6,7). 

We take the best values for Am32 known at present 
(see for example Ref.[ll] and references therein) and from 
CHOOZ evidence [34, 35] the value of the (13) angle: 1.3 x 

10-3 eV^ < \AmltJ < 3.1 x 10-3eV^ 0.90 < sin^ 2^23 < 
1.0, sin^ 013 < 0.047, 90% C.L.. With this values and for 
those values obtained previously for the 2x2 Mq matrix 
(Eq.[6,7]) we finally obtain an estimation for the three 
neutrino squared mass matrix (10-^ eV^ units) : 



8.1 ±6.5 
8.4 ±6.6 
-27.7 ±28.0 



8.4 ±6.6 
8.6 ±6.9 
-27.6 ± 29.0 



-27.7 ±28.0' 
-27.5 ± 29.0 
202.0 ± 60.0 , 



One can go further supposing a concrete value for the 

free parameter from elsewhere. If wc take m2 >> Am^ 
then we can directly write the mass matrix 



M = m2l + 



2m2 



(9) 



Supposing for example mi 



M 



1 + (4.0 ± 3.2) 10" 
4.2 ± 3.2 lO-'^ 
13.5 ± 14.0 lO-'^ 



1 eV, we obtain {eV units) 

-13.5 ±14.0 10" 



4.2 ± 3.2 10-" 
1 ± (4.3 ± 3.5) 10 
-13.5 ± 14.5 10" 



-13.5 ± 14.5 10 



1 



(100.0 ±30.0) 10" 



this is now the final three neutrino mass matrix which 
can be obtained from present oscillation evidence coming 
from solar and reactor neutrinos. 



5 Summary and Conclusions 

Wc have presented an up-to-date analysis including the re- 
cent KamLAND results, the SNO-phase II spectrum and 
all other solar neutrino data The active neutrino oscilla- 
tions hypothesis has been confirmed, and the decoupling of 
the atmospheric zim^-solar Am^ justifies a 2-fiavour anal- 
ysis as the one presented here. This justification is even 
stronger if we have into account the large experimental 
disparity among solar, earth reactor and atmospheric evi- 
dence and the very much different accuracy which can be 
obtained in each of them for the parameters of the /xr and 
e/Lt neutrino sectors. Moreover, the consideration in the 
analysis of the atmospheric data would only slightly mod- 
ify the best values and allowed regions for the parameters. 
These modifications would be well within the error bars of 
these paramters according to the present determination. 
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The results presented along this work show how due 
to the increased statistics, the inclusion of the new Kam- 
LAND data determines with good accuracy the value of 
AniQ, clearly selecting the LMAI solution, and brings us 
to a new era of precision measurements in the solar neu- 
trino parameter space [32]. 

Wc have introduced in this work diverse novelties in 
the treatment of the SNO and Kl spectra. For the first 
one we have improved upon previous works in the full 
consideration of the sistematic correlations. For the KL 
spectrum we have studied the variability of the best fit re- 
sults with respect the statistical method in use. We have 
shown that appreciable differences can be obtained. We 
believe that a careful study and proper statistical treat- 
ment of the KL evidence is needed. Significant differences 
on the values of the oscillation parameters can be obtained 
basically due to poor statistics. These apparition of these 
differences can be easily missed or obscured by analysis 
which include large quantities of diverse data without the 
needed care of the individual components. 

We have obtained the allowed area in parameter space 
and individual values for Z\m^ and tan^ with error esti- 
mation from the analysis of marginal likelihoods. We have 
shown that it is already possible to determine at present 
active two neutrino oscillation parameters with relatively 
good accuracy. In the framework of two active neutrino 
oscillations we obtain 

Am^ = 8.20 ± 0.08 x lO'^eY^, tan^ 6 ^ 0.50ta l,l. 

We estimate the individual elements of the two neutrino 
mass matrix, we show that individual elements of this ma- 
trix can be determined with an error ~ 10% from present 
experimental evidence. 

The use of the SNO phase-II spectrum in the data set 
has mainly two effects: 1) a slight reduction in the overall 
area of the exclusion plot and 2) a slight decrease in the 
best-fit Aw?Q. 

The decrease in the best-fit mass squared difference 
can be understood by the fact that by including the SNO 
spectrum, we increase the statistical relevance of solar 
neutrino data, which prefer smaller Am?. Furthermore, 
the oscillation pattern (whose information is contained in 
the spectrum) is more sensitive to Arn^. 

It is interesting to note that the KamLAND data alone 
still continue to predict, for both their analyses, a value 
of tan^ 9 smaller than the one obtained with the previous 
data, and significantly different from 1, consequently mak- 
ing the aesthetically pleasing bi-maximal-mixing models 
strongly disfavored. This result confirms what was already 
evident in the solar neutrino data analyses. Nevertheless, 
improvement on the determination of tan^ 9 is necessary 
and it is known that KamLAND is only slightly sensi- 
tive to this mixing parameter. The (lower) accuracy with 
which we determine the solar mixing angle is evident in 
the marginalized likelihood plots of fig. [2]. Planning of 
future super-beam experiments aimed at determining the 
^13 and eventual CP violating phases relies on the most 
accurate estimation of all the mixing parameters [33] . It is 
expected that future solar neutrino experiments, notably 



phase-II SNO (higher statistics, due to be made public 
soon) and eventually future low energy experiments, and 
phase-Ill SNO (with helium) will further restrict the al- 
lowed range of parameters. 
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Table 1. Summary of Kamland spectrum information extracted from [1]. Relative statystical errors only are reported. 



T.fs (MeV) 


Evnts./500 keV 


Tefs (McV) 


Evnts./500 keV 


5.5- 6.0 


225 


9.5-10.0 


155 


6.0- 6.5 


225 


10.0-10.5 


80 


6.5- 7.0 


220 


10.5-11.0 


95 


7.0- 7.5 


255 


11.0-11.5 


55 


7.5- 8.0 


235 


11.5-12.0 


40 


8.0- 8.5 


225 


12.5-13.0 


15 


8.5- 9.0 


155 


13.5-14.0 


15 


9.0- 9.5 


145 


14.5-15.0 


5 


T.fs (MeV) 


Evnts./500 keV 


T.fs (MeV) 


Evnts./500 keV 


5.5- 6.0 


840 ± 20 


9.5-10.0 


180 ± 10 


6.0- 6.5 


785 ± 20 


10.0-10.5 


110 ±10 


6.5- 7.0 


705 ± 20 


10.5-11.0 


110 ± 10 


7.0- 7.5 


680 ± 20 


11.0-11.5 


50 ±10 


7.5- 8.0 


560 ± 15 


11.5-12.0 


40 ± 10 


8.0- 8.5 


470 ± 15 


12.5-13.0 


10± 10 


8.5- 9.0 


245 ± 10 


13.5-14.0 


10 ±5 


9.0- 9.5 


205 ± 10 


14.5-15.0 


0±5 



Table 2. (top) CC Energy spectrum observed at SNO (taken from fig. 36 of ref [8]). (bottom) SNO neutrino spectrum ( fig.24 
from Ref. [8]. 



Experiment 


References 
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SNO 
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CHOOZ 


[34] 




Palo Verde 


[35] 




KamLAND 
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Table 3. References from where we draw 


the data used in 


our analysis. 






tan^6l 


from X'' minimization 






KL (Sp-l-Gl)-FSolar -|- SNO (Sp) 


7.89 X 10"^ 


0.40 


KL (Sp-l-GlHSolar -|- SNO (Gl) 


8.17 X 10"^ 


0.40 


from marginalization 






KL (Sp-l-Gl)-FSolar -|- SNO (Sp) 


8.2in X 10-^ 


U.OU_g Qg 



Table 4. Mixing parameters from minimization and likelihood marginalization. 
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Fig. 1. (left) Allowed region in the (tan^ 6, ArriQ) plane for the global analysis, which includes the previous solar data (see e.g. 
[6] for details) and all KamLAND Global results, (right) Best fit solution for the spectrum analysis, including all previous solar 
data, short baseline reactor data and the KamLAND spectrum. Best fit point given in table [4]. 




Fig. 2. Marginalized likelihood distributions for each of the oscillation parameters Attiq (left) and tan^ 6 (right) corresponding 
to the totality of solar and KamLAND data. The curves are in arbitrary units with normalization to the maximum height. 



